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Abstract  In  this  paper,  we  concentrate  on  one  of  the 

The  paper  presents  a  study  of  stability  of  most  interesting  features  of  explosive  boiling, 

explosive  boiling,  where,  in  addition  to  surface  which  is  a  very  high  evaporation  rate  attainable  in 

evaporation,  a  vapor  bubble  grows  within  a  this  process.  One  cause  for  such  strong 

highly  superheated  liquid  droplet,  immersed  in  a  evaporation  is  thought  to  be  the  significant 

liquid  or  gas  medium.  To  get  better  insight  into  increase  of  the  area  of  the  evaporation  surface 

the  problem  two  simpler  but  related  problems  are  caused  by  instability  of  the  bubble  interface2'3,5, 

studied  before  the  full  stability  problem  is  This  instability,  first  observed  by  Shepherd  & 

treated.  First,  the  stability  of  an  evaporating  Sturtevant3,  manifests  itself  as  wrinkling  and 

highly  superheated  liquid  droplet  is  analyzed,  in  roughening  of  the  vapor  bubble  surface  followed 

order  to  estimate  the  influence  of  the  outer  by  its  distortion.  The  increased  area  of  the 

evaporation  from  the  droplet  surface.  The  linear  evaporation  surface  provides  the  necessary  heat 

stability  of  the  process  at  the  final  stages  of  transfer  to  support  explosive  boiling, 

explosive  boiling,  when  the  droplet  forms  an  Recently  Shusser  &  Weihs6  proposed  a 

expanding  liquid  shell,  is  studied  next.  Finally,  mathematical  model  describing  growth  of  an 

the  general  case  of  explosive  boiling  stability  is  internal  vapor  bubble  produced  by  homogeneous 

considered.  It  is  shown  that  the  process  is  nucleation  within  a  liquid  droplet  during 

unstable  as  indeed  has  been  found  in  existing  explosive  boiling.  The  predictions  of  the  model 

experiments.  were  confirmed  by  existing  experimental  results 

for  explosive  boiling  of  superheated  droplets3,5. 

Introduction  The  instability  of  explosive  vapor  bubble  growth 

The  main  task  of  an  aircraft  engine  is  to  was  n°f  however,  considered  by  Shusser  & 

transform  the  chemical  energy  of  the  fuel  into  the  Weihs6.  Proper  understanding  of  this  instability 

kinetic  energy  of  the  air.  In  air-breathing  jet  W^1  throw  light  on  the  physics  of  explosive 

engines,  this  is  achieved  by  combustion  of  liquid  boiling. 

fuel  in  a  combustion  chamber.  For  this  chemical  Explosive  boiling  instability  is  related  to  the 

reaction  to  occur,  the  liquid  fuel  must  be  injected  instability  of  laminar  flames  discovered  by 

into  the  airstream,  atomized,  and  the  vapor  must  Landau7  and  investigated  for  spherical  flames  by 

be  mixed  with  the  air.  Studying  the  behavior  of  Istratov  &  Librovich  .  It  is  also  connected  to  the 

fuel  droplets  injected  into  a  hot  airstream  is  instability  of  evaporation  surfaces  and  to  the 

important  for  proper  designing  of  combustors  for  problem  of  spherical  bubble  stability  .  The 

air  breathing  jet  propulsion  engines1.  unperturbed  state  is  time-dependent  (growth  of  a 

To  achieve  high  efficiency,  the  jet  engine  spherical  vapor  bubble  within  a  liquid  droplet) 

designer  generally  attempts  to  maximize  the  air  and  hence  normal-mode  analysis  is  not 

flow  rate  and  the  temperature  in  the  combustion  appropriate. 

zone.  When  immersed  in  a  hot  combustion  gas,  Instabilities  developing  on  the  outer  surface 

liquid  fuel  droplets  are  suddenly  and  drastically  liquid  droplet  can  be  capillary  instability  if 

heated.  For  high  gas  temperatures  and  high  flow  tbe  droplet  is  situated  in  a  liquid  medium  or 

rates,  as  in  modem  aircraft  engines,  this  may  evaporating  surface  instability  if  the  medium  is 

cause  high  superheating  of  the  droplets.  gaseous.  To  analyze  the  interaction  of  the  process 

The  process  of  rapid  phase  transition  from  on  both  interfaces  one  must  consider  two  types  of 

highly  superheated  liquid  to  vapor  is  called  perturbations  which  we  shall  call  symmetric 

explosive  boiling2'4.  It  is  usually  so  fast  and  and  "antisymmetric"  in  analogy  from  stability  of 

violent  that  it  resembles  an  explosion.  Explosive  liquid  films16  or  annular  liquid  jets  7. 

boiling  of  fuel  droplets  inside  the  combustion 
chamber  can  damage  the  engine. 
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The  physical  mechanism  responsible  for  this 
instability  is  not  fully  understood  at  present. 
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Sturtevant  &  Shepherd18  used  the  Landau  theory 
to  estimate  stability  limits  and  growth  rates  for 
explosive  boiling.  Recently  Lee  &  Merte19 
showed  that  approximating  the  evaporation 
surface  as  plane  and  using  instantaneous 
information  for  a  growing  spherical  vapor  bubble 
one  can  reasonably  predict  the  occurrence  of 
instability  and  its  wavelength.  Nevertheless,  the 
understanding  of  rapid  evaporation  instability  for 
the  spherical  case  has  not  been  achieved  yet19. 

Our  aim  is  to  investigate  the  linear  stability 
of  explosive  boiling  of  a  liquid  droplet  in  liquid 
or  gas  medium.  We  concentrate  on  hydrodynamic 
aspects  of  the  problem  but  consider  spherical 
geometry  and  include  both  interfaces  into  the 
analysis.  We  start  by  studying  two  simpler  but 
related  problems  before  treating  the  general  case. 

The  plan  of  the  present  paper  is  as  follows. 
First,  the  stability  of  evaporation  of  a  highly 
superheated  liquid  droplet  is  studied.  Next,  the 
problem  of  linear  stability  of  a  thin  expanding 
liquid  spherical  shell  is  considered.  Then  we 
proceed  to  analyze  the  general  case  of  explosive 
boiling  stability. 


R0  =  R(0)  and  p,pv  are  the  densities  of  the 
liquid  and  the  vapor,  respectively; 
a  =  pv  /  p  (0  <  a  <  1).  Furthermore  a  is  the 
surface  tension;  r  is  the  radial  coordinate  in  the 
spherical  coordinates  (r,0,(p)  ;  t  is  the  time. 

For  spherical  incompressible  flow  with 
constant  evaporation  rate  to  satisfy  the 
conservation  of  mass,  the  droplet  radius  must 
decrease  linearly  with  time 

J 

R  =  R0— t  (i) 

P 

From  the  conservation  of  mass  and 
momentum  at  the  interface,  one  can  calculate  the 
pressure  and  velocity  fields  (see  Shusser21). 

The  droplet  is  now  assumed  to  have  the 
following  surface  shape 

r  =  R(t)  +  e(0;<p;t)  (2) 

where  8  «  R  (see  insert  in  Fig.  1). 

We  denote  the  components  of  perturbation 
velocity  in  spherical  coordinates  and  perturbation 

of  pressure  Vf_ ,  VQ_ ,  ,p_  on  the  droplet 


Stability  of  a  Highly  Superheated  Evaporating 
Liquid  Droplet 

Take  a  spherical  highly  superheated  liquid 
droplet  surrounded  by  vapor  of  the  same 
composition.  The  droplet  evaporates  creating 
vapor  flow  in  the  host  medium  as  shown  in  Fig. 
1.  Here  no  internal  vapor  bubbles  are  produced. 
Our  purpose  is  to  study  the  linear  stability  of  this 
process. 

The  stability  of  evaporation  has  been 
investigated  only  for  evaporation  from  plane 
surfaces9'11,20.  Moreover,  most  of  the  work  has 
been  devoted  to  studying  marginal  stability, 
which  may  be  inappropriate  for  investigation  of 
the  stability  of  evaporation  (Prosperetti  & 
Plesset11,  p.  1590).  Therefore,  as  a  first  stage  in 
the  study  of  the  stability  of  an  explosively 
evaporating  droplet  we  make  several  simplifying 
assumptions  based  on  the  observations  of  high 
superheating  and  very  strong  evaporation.  These 
assumptions  are  consistent  with  explosive 
boiling. 

We  assume  that  the  vapor  and  the  liquid  are 
inviscid  incompressible  fluids,  that  evaporation 
rate  in  the  base  flow  is  constant,  that  the  flow 
perturbations  do  not  influence  the  evaporation 
rate,  and  that  the  flow  field  is  spherically 
symmetric. 

Let  U+,p+  andu_,p_  be  the  velocity  and 
the  pressure  on  the  vapor  and  liquid  phase, 
respectively.  p0  is  the  pressure  far  from  the 
droplet,  J  is  the  evaporation  rate  (p.u.  surface 
area),  R=R(t)  is  the  droplet  radius  where 


and  Vf+ ,  V0+  ,p+  in  the  vapor.  At  the 

interface  they  satisfy  boundary  conditions  which 
are  the  conservation  of  mass,  the  constancy  of  the 
evaporation  rate  and  the  conservation  of  three 
components  of  momentum. 

Istratov  &  Librovich8  solved  a  similar  (but 
inverse)  problem  of  spherical  flame  stability. 
Following  their  solution  we  obtain  that  at  the 
unperturbed  surface  of  the  droplet  r  =  R(t) 
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dr  -  8t 
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where  A  is  the  perturbed  droplet  surface 
curvature. 

For  inviscid  fluid,  the  flow  inside  the 
perturbed  droplet  remains  irrotational.  Defining 

the  perturbation  potential  0_  and  choosing  an 
appropriate  solution  of  the  Laplace  equation 

<J>:=f(t)r”Ynm(0;(p)  (8) 


(3) 

(4) 

(5) 

(6) 
(7) 
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where  Ynm  are  spherical  harmonics  one  obtains 
the  components  of  perturbation  velocity  and 
perturbation  pressure 
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df  „  m 

p-=-pdFr  Y" 

(10) 

The  flow  in  the  vapor  phase  will  be 
rotational  due  to  the  creation  of  vorticity  at  the 
distorted  surface  of  the  perturbed  droplet.  It 
should  be  stressed  that  it  is  not  possible  to  use  the 
Laplace  equation  for  pressure  perturbations  in  the 
vapor  phase7’8  because  there  is  a  vapor  flow  in 
the  unperturbed  case.  It  turns  out  that  our 
problem  (outer  evaporation  from  the  droplet 
surface)  is  much  more  difficult  than  the  problem 
of  Istratov  &  Librovich8  (expansion  of  spherical 
flame)  or  the  problem  of  evaporation  into  an 
inner  bubble. 


Due  to  the  velocity  perturbation  V+  being 
solenoidal,  we  can  divide  it  into  toroidal  and 
poloidal  parts  T,S  (Chandrasekhar22,  p.  225). 
The  vorticity  £2  can  also  be  divided 


£2  =  T  +  S  (11) 

where  the  functions  T,  S  satisfy  (Ref.  22,  p. 
226) 

~  n(n  +  l)  d2S 

T  = - ^ - S-TT  (12) 

r  dr 

S  =  T  (13) 

The  radial  component  of  vorticity  is 

n(n  +  l) 

Qr  =  — — ^  ~T(r)Yn  (14) 

We  now  show  that  in  the  linear  approximation 
£2r  vanishes. 

Indeed,  linearizing  the  vorticity  equation  for 
inviscid  flow  we  obtain  for  £2r 


9£2 


3£2r 


£2„ 


dir 


■  =  0  (15) 


3t  +  9r  dr 

and  therefore  £2r  is  not  identically  zero  only  if  it 
is  created  at  the  perturbed  interface  (2).  On  the 
other  hand,  the  vorticity  is  of  order  O(s)  and 


therefore  neglecting  second-order  terms  one  can 
calculate  it  at  the  surface  of  the  unperturbed 
droplet  r  =  R .  Then  using  the  definition  of 
vorticity  and  boundary  conditions  (5-6)  one 
obtains  that  the  radial  component  of  vorticity  is 


continuous  at  the  interface.  It  is  zero  inside  the 
droplet,  so  that  it  must  vanish  in  the  vapor  too. 

We  have  therefore  obtained  that  £2r  is  a 
second-order  quantity  and  hence  negligible  in  the 
investigation  of  linear  stability.  Thus  we  put 

T(r)  =  0  (16) 

That  is,  in  the  vapor  phase  the  flow-field  has  only 
a  poloidal  part  and  the  vorticity  field  has  only  a 
toroidal  part.  This  is  in  accordance  with 
Prosperetti  (Ref.  23,  p.  344)  who  stated  that  in 
perturbed  spherical  flows  a  poloidal  part  of  the 
vorticity  cannot  be  generated  if  it  vanishes  at  the 
initial  moment. 

Returning  to  the  calculation  of  the  velocity 
field  within  the  vapor,  relation  (12)  is  now  an 
ordinary  differential  equation,  with  general 
solution 

C(t) 

S=-^Z  +  D(t)r"-1+Sp  (17) 

where  Sp  is  the  particular  solution. 

We  assume  that  when  r  — »  the  vorticity 
tends  to  zero  sufficiently  rapidly  so  that 
S  — »  0  too.  Then 

D(t)  3  0  (18) 

We  are  interested  in  the  behavior  of  S(r;t) 
near  the  droplet  surface.  Therefore  we 
approximate  T(r;t)  by  its  value  at  the  droplet 
surface,  which  we  denote  F(t) 

T(r;t)  ~  T(R(t);t)  =  F(t)  (19) 


Then  calculating  Sp  we  obtain 

_  C(t)  F(t)r2 
rn  +  n(n  +  l)-2 

The  solution  (20)  always  exists,  as  we  are 
interested  in  the  perturbation  modes  for  which 
n  >  2 . 

From  (16)  and  (20)  one  obtains  the  perturbed 
velocity  field  in  the  vapor,  while  the  pressure 
perturbation  can  be  calculated  from  the  linearized 
Euler  equation21.  In  addition,  the  droplet  surface 
distortion  is 


(20) 


8  =  a,(t)Ynm(0;(p)  ,  n  >  2  (21) 

Substituting  the  pressure  and  velocity  fields 
into  the  boundary  conditions  (3-7),  we  can  obtain 
the  equation  for  the  amplitude  of  the  droplet 
surface  perturbation  a^t)  (see  Shusser21).  It 
follows  from  the  solution  of  this  equation  that 
when  R  — >  0 

a  2fci 

^~R  4  (22) 

Here 
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(n-1) 
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3  — 


n  +1 


-  +  a 


(23) 


One  sees  that  P  >8/3  when 
n  >  2  and  0  <  a  <  1  and  therefore 


aj  /  R  — »  0  when  R  — >  0 ,  i.e.,  droplet 
surface  perturbations  tend  to  zero  faster  than  the 
droplet  radius.  The  solution  is  thus  stable. 

A  question  arises  about  the  validity  of  this 


analysis,  due  to  the  assumption  of  “  «  1 


(equation  2)  which  is  only  justified  a  posteriori 
here. 

However,  returning  to  our  main  objective, 
stability  of  explosive  boiling,  the  results  of  this 
section  enable  us  to  ignore  the  outer  evaporation 
from  the  droplet  surface  when  analyzing  the 
stability  of  a  liquid  droplet  boiling  explosively,  as 
R  never  tends  to  zero  in  the  full  problem,  where 
an  internal  vapor  bubble  is  produced. 

The  decay  is  faster  for  higher  wave  numbers 
and  when  vapor  and  liquid  densities  are  closer  in 
value  (a  — »  1) ,  as  seen  from  (23). 


It  is  interesting  to  compare  our  results  for 
evaporating  droplets  and  the  results  for 
expanding  spherical  flames8  with  those  for 
spherical  bubbles12'15.  In  the  latter  case  there  is 
no  flow  through  the  bubble  surface  and  therefore 
the  perturbed  flow  remains  irrotational. 

One  sees  from  the  comparison  that  the 
results  are  opposite.  The  growing  bubble 
(R  —>  oo)  is  stable  and  the  collapsing  bubble 
( R  — >  0 )  is  unstable.  On  the  other  hand,  the 
expanding  flame  (R  — »  oo)  is  unstable  and 
evaporating  droplets  ( R  — »  0 ),  as  we  have  just 
shown,  are  stable. 


and  the  inner  bubble  pj  is  constant  and 
uniform. 

The  last  assumption  corresponds  to  the 
existence  of  very  weak  evaporation  into  the  inner 
bubble.  Because  it  is  weak  one  can  neglect  its 
influence  on  the  shell  mass  and  the  stability,  but 
due  to  the  low  density  of  the  vapor  it  is  sufficient 
to  support  constant  pressure  within  the  bubble. 

We  utilize  the  assumption  of  a  thin  shell  by 

neglecting  quantities  of  order  h2  /R2  both  in 
the  base  solution  and  in  the  perturbed  flow. 

Let  4ttM  be  the  mass  of  the  shell.  Then 
from  the  conservation  of  mass 


4ttM  =  4nhR2p 


1  +  - 


12R2 


and  after  neglecting  the  second-order  quantities 

2  M 

hR  =  —  =  const 
P 

To  this  approximation  the  flow-field  in  the  shell 
is  (the  dot  denotes  a  time  derivative) 

R2 


v  =  * 


-R 


Then  from  the  unsteady  Bernoulli  equation  of 
hydrodynamics24  and  the  boundary  conditions  at 
both  surfaces  one  obtains  in  the  linear 
approximation  the  equation  for  the  shell  radius 


4°  (Pj-pJ 


r+TTr~ 

M 
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R2  =0 


(24) 


(25) 


(26) 


(27) 


This  equation  does  not  include  a  R  term.  This 
means  that  there  is  no  damping  term  in  the 
equation  and  therefore  thin  shell  oscillations  do 
not  decay  or  amplify  in  the  linear  approximation. 
Integrating  (27),  we  obtain  for  the  expanding 

shell 


Stability  of  a  Thin  Expanding  Spherical  Shell 
Consider  a  thin  spherical  liquid  shell  of 
density  p  and  surface  tension  a  expanding  in  a 
gas  medium  as  shown  in  Fig.  2.  The  shell  is 
characterized  by  its  mean  geometric  radius  R(t) 
and  its  thickness  h(t).  This  model  describes  the 
late  stage  of  explosive  boiling. 

The  instability  of  an  expanding  liquid  shell  is 
stronger  for  boiling  in  a  gas.  This  results  from  the 
large  difference  in  density  between  liquid  and 
gas,  which  facilitates  the  motions  of  shell 
segments.  Therefore  we  analyze  this  case. 

We  assume  that  the  shell  is  thin  relative  to  its 
radius  (h«R),  that  the  unperturbed  (base) 
solution  is  spherically  symmetric,  that  the  liquid 
is  inviscid  and  incompressible,  that  gas  and 
internal  vapor  density  are  negligible,  that 
evaporation  from  the  outer  surface  of  the  droplet 
is  negligible,  and  that  the  pressure  in  the  host  gas 


where 


R  =  J-bR3  -aR2  +c 


4o  ,  Pi-P* 


and  c  is  defined  by  initial  conditions. 

The  expression  within  the  radical  in  (28) 
must  be  positive.  This  can  be  considered  as  a 
condition  on  the  shell  radius  R  for  sustaining  a 
thin  shell. 

We  perturb  the  base  solution  so  the  shell 
radius  and  the  shell  thickness  are  now 


R(t)  +  r|(t;0;(p)  and  h(t)  +  £(t;0;(p) . 


The  perturbations  are  small  and  therefore 
T|  «  R  and  «  h  and  the  terms  of  order 


0(n2)and  0(£2  )  are  negligible. 


(28) 


(29) 


4 


The  perturbed  shell  is  situated  between  two 
slightly  perturbed  spheres.  Their  equations  are 

h  5 

r  =  R  +  -  +  -n  +  -  (30) 

From  conservation  of  mass 

4ftM  =  p(v2  -  V,)  (31) 


where  Vj  and  V2  are  the  volumes  of  the  inner 
and  outer  perturbed  spheres,  respectively.  To 
linear  approximation 


Substituting  (25),  (32-33)  into  (31)  and 

neglecting  the  second-order  terms  in  h  we  obtain 


f2,t  f  (R2£  +  2Rhr|)  d0dcp  =  0  (34) 

Jo  Jo 

One  can  conclude  that  there  exist  two  types  of 
perturbations  (see  Fig.  2).  If  the  perturbation  of 
the  radius  r\  is  not  identically  zero  then 
t  2h 

(35) 

The  perturbations  of  the  shell  radius  are 
analogous  to  the  antisymmetric  (sinuous) 
perturbations  observed  in  plane  liquid  films16  or 
annular  liquid  jets17.  For  spherical  flows,  they 
also  cause  perturbations  of  the  shell  thickness, 
given  by  (35). 

It  is  possible  to  perturb  the  thickness  leaving 
the  shell  radius  unchanged  provided  the  shell 
volume  does  not  change,  i.e.,  if 

J'2n  rn 

0  J0  k  dedtp  =  0  (36) 

This  type  of  perturbation  is  analogous  to  the 
symmetric  (varicose)  perturbations  for  the  plane 
film16.  Any  perturbation  can  be  represented  as  a 
linear  combination  of  symmetric  and 
antisymmetric  ones.  Therefore  each  type  of 
perturbation  can  be  analyzed  separately.  We 
begin  the  analysis  for  the  antisymmetric  case. 

The  perturbation  potential  is  a  solution  of  the 
Laplace  equation  and  therefore  the  full  potential 
of  the  perturbed  flow  O  is  given  by 

O  =  -^  +  f^  +  a2(t)rnVm(0;cp)  (37) 


Assuming  for  the  perturbation  T|  a  solution  of 
the  form 

Tl  =  a3(t)Ynm(e;cp)  (38) 

we  can  find  the  functions  a1(t),a2(t),a3(t) 

from  the  boundary  conditions,  which  are  two 
kinematic  boundary  conditions  (Lamb24,  p.  7)  and 
the  conservation  of  momentum21. 

We  are  interested  in  the  behavior  of  the 
solution  for  t  — >  °°,  i.e.  R  — >  00 .  Neglecting 

the  terms  of  order  R  ~2 ,  one  obtains  (see 
Shusser21)  that  when  R  00 

a.  Ci  in 

Y~^t+c2r  <39) 

For  stability  a3  /  R  must  remain  bounded. 
Therefore  one  of  the  solutions  is  always  unstable 
and  the  conclusion  is  that  the  expansion  of  a  thin 
liquid  shell  is  unstable  for  antisymmetric 
perturbations.  It  can  be  noted  that  the  instability 

is  relatively  weak,  growing  as  R 1/2 . 

Considering  symmetric  perturbations,  we 
first  mention  that  they  are  perturbations  of  the 
shell  thickness.  They  are  initially  much  smaller 
than  the  thickness  itself  even  in  the  thin-shell 
approximation.  Therefore,  in  analyzing  stability 
of  symmetric  perturbations  we  retain  the  terms  of 

order  0(h2  /  R2  )  in  the  base  solution. 

Repeating  the  previous  calculations 
including  the  0(h2  /  R2)  terms,  we  obtain  the 
shell  thickness  h  and  the  flow  potential  O  are 
now 


M2  | 

12p2R6J 

(40) 

3h2  ^ 

1_  4R2; 

(41) 

Proceeding  as  previously,  we  look  for  a 
solution  for  the  thickness  perturbation  £  of  the 
form 


£  =  a3(t)Ynm(e;(p) 

with  the  potential  in  the  perturbed  flow 

.  R2  J.  3h2  ) 

4>  = - R  1 - -  + 


a2(t)rn 


(42) 


(43) 


Then  the  equation  for  a3(t)  is 

d2a,  da3 

d#  +  P(R)dR  +Q(R)a3=°  (44> 

where  P(R)  and  Q(R)  are  given  by  Shusser21. 

From  the  leading  terms  of  P(R),  Q(R)  for  R 
sufficiently  large 
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Q  >  0  ;  —  +  2PQ<0  (45) 

dK 

and  therefore  (Birkhoff13)  each  solution  of  (44) 
tends  to  infinity  when  R  — >  °c> . 

To  find  the  behavior  of  the  solution  for 
R  — >  00 ,  we  substitute  in  Eq.  (44)  the  first 
terms  in  the  asymptotic  expansions  of  P  and  Q 
for  large  R.  Comparing  the  behavior  of  the  shell 
thickness  h  and  the  amplitude  of  its  perturbation 
a  3 ,  one  sees  that  when  R  — >  00  and  h  — »  0 

h*13'8  ->°o  (46) 

h  h->0 

That  is,  there  is  strong  instability  for  all  wave 
numbers. 

We  obtained  that  both  perturbations  of  the 
shell  radius  (antisymmetric)  and  perturbations  of 
the  shell  thickness  (symmetric)  are  unstable. 
When  the  shell  radius  R  tends  to  infinity  and  the 
shell  thickness  h  tends  to  zero,  the  appropriate 
perturbations  normalized  by  R  or  h,  respectively, 

are  asymptotic  to  VR  and  h~13/8 . 

Squire16  has  shown  that  for  thin  film  stability 
both  types  of  perturbations  are  unstable.  In 
variance  with  our  solution,  however,  in  the  plane 
case  antisymmetric  perturbations  grow  much 
faster  than  symmetric  ones.  For  a  spherical  shell 
we  obtained  weak  instability  for  perturbations  of 
the  shell  radius  and  strong  instability  for  those  of 
the  thickness. 

One  sees  that  the  fact  that  both  shell  radius 
and  shell  thickness  are  time-dependent  in  our 
spherical  case  does  not  change  the  general 
conclusion  about  stability  but  changes  its  nature 
because  the  most  dangerous  perturbations  are 
now  of  different  type.  The  weak  instability 
obtained  for  the  antisymmetric  perturbations  may 
tie  in  with  the  related  process  of  gas  bubble 
growth  with  a  constant  pressure  inside  the 
bubble,  which  is  stable12"15. 

The  goal  of  this  section  was  to  identify  the 
most  dangerous  perturbations  for  stability  of 
explosive  boiling.  These  are  the  symmetric 
(varicose)  perturbations. 

Stability  of  Explosive  Boiling 

The  physical  situation  is  depicted  in  Fig.  3. 
A  vapor  bubble  of  radius  R^t)  grows  within  a 
highly  superheated  liquid  droplet  of  radius 
R2(t) ,  which  is  situated  in  liquid  or  gas 
medium.  For  simplicity,  we  assume  that  the 
bubble  is  situated  in  the  center  of  the  droplet. 
Though  an  approximation,  this  assumption  was 
shown  reasonable  for  a  broad  range  of  physical 
situations6. 

We  assume  that  all  the  fluids  are  inviscid  and 
incompressible  and  the  evaporation  rate  is 


constant  and  is  not  influenced  by  the 
perturbations.  We  also  neglect  the  flow  inside  the 
vapor  bubble.  Utilizing  the  results  of  the  previous 
sections  we  neglect  the  effects  of  possible 
evaporation  from  the  outer  surface  of  the  droplet 
and  consider  only  symmetric  perturbations. 

Due  to  the  constant  rate  of  evaporation,  the 
rate  of  growth  of  the  bubble  radius  will  also  be 
constant3.  Denoting  this  constant  rate  U  we  write 


=Ut 

and  using  conservation  of  mass 


R2 

R0 


/ 

l  +  (l-a) 

v 


Ur 

R0> 


(47) 


(48) 


Here  R0  is  the  initial  radius  of  the  droplet; 
a  =  pv  /  p;  p,pv  are  the  densities  of  the 
droplet  liquid  and  the  vapor,  respectively. 

We  add  small  symmetric  perturbations 
8(0;<p;t)(8  «  Rj,R2)  to  both  interfaces  so 
their  equations  are  now 

ii  =  Rx  (t)  +  8(0;cp;t)  (49) 

r2  =R2(t)-e(0;<p;t)  (50) 

From  conservation  of  mass 

r 2n  fn 

Jo  Jq  8(0;(p;t)sin0d0d(p  =  0  (51) 

Distortion  of  the  droplet  shape  causes 
perturbations  of  the  velocity  Vr.,Ve.  ,V9i  anc* 


the  pressure  Pj  (i=l,2,3  for  the  vapor  bubble,  the 
liquid  droplet  and  the  host  fluid,  respectively). 
These  perturbations  satisfy  the  boundary 
conditions  at  the  surface  of  the  bubble  and  the 
surface  of  the  droplet. 

We  write  the  conditions  of  conservation  of 
mass,  conservation  of  the  tangential  components 
of  momentum  and  constancy  of  the  evaporation 
rate  at  the  surface  of  the  unperturbed  bubble 
r  =  Rj  as 


v  r  =-^T 

ri  at 

(52) 

.  avrj  3e 

vr  +8  _  -  - 

2  3r  3t 

(53) 

1  3s 

v„  =  v„  +  — —  V, 

1  2  R,  30  2 

(54) 

1  3s 

y  =  v  H - V 

*2  R,sin0  3(p  f2 

(55) 

The  kinematic  conditions  at  the  surface  of  the 

unperturbed  droplet  r  =  R2  are 

,  3e  3V 

V.  +  -  6  . 

f2  3t  3r 

(56) 
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8s  8V 


V_  +—  =  8- 
r’  3t  3r 

Conservation  of  the  normal 

momentum  at  r  =  R }  and  at  r 

the  dynamic  boundary  condition 


“Pi 


r=R, 


(57) 

component  of 
=  R2  results  in 


V  R'  Rv 


(  ■  ;)p.  1 

.  3P,) 

HirJ 

r=R2 

Lp=+e^rJ 

(58) 


r=R, 


Here  Vr2 ,  Vr3 ,  P2 ,  P3  denote  the  unperturbed 


solution. 

We  write  the  perturbation  e  as,  using  (51) 
e(0;(p;t)  =  f4(t)Y™(0;cp)  (n*0)  (59) 


The  perturbed  pressure  and  velocity  fields  in 
the  vapor  bubble,  the  liquid  droplet  and  the  host 
liquid  were  calculated  by  Shusser21.  By 

substituting  them  into  the  boundary  conditions 
(52-58)  the  following  equation  for  the  amplitude 
of  the  shape  perturbation  is  obtained 

d2f,  df. 

A(t)— f  +  B(t)-f  +  C(t)f4=0  (60) 
dt  dt 

The  functions  A(t),  B(t),  C(t)  are  written  out 
in  Ref.  21. 

We  start  the  analysis  by  defining  the 
conditions  for  stability.  The  process  of  explosive 
boiling  occurs  during  a  finite  time  until  the 

droplet  evaporates  completely,  i.e.,  until 


Rn 


1 "  tf  "  Ua1'3 


(61) 


We  propose  to  analyze  the  stability  by 
investigating  the  behavior  of  the  solution  of  (60) 
when  t— >tf.  It  is  possible  that  the 
perturbations  grow  so  fast  that  the  droplet  is 
destroyed  before  t  f .  Therefore  our  method  can 
prove  only  instability  of  the  process. 

Taking  new  variables 


tf  -  t  ;  X  =  ■ 


Ua 


1/3 


a 


1/3 


R„ 


y  =  T"f4  (62) 


one  can  write  the  approximate  form  of  (60)  in  the 
limit  T— >0  as 

d2y  dy 

x(a,x  +  a2)  — f  +(b,x  +  b2)  —  + 
dx  dx 


(63) 


(c,x  +  c2)y  =  0 

The  coefficients  ai,  a2,  bh  b2,  Ci,  c2  are  given  by 
Shusser21. 


No  general  analytic  solutions  for  Eq.  (60)  are 
known.  Studying  the  behavior  of  the  solution 
when  X  — >  0  by  using  the  method  of 
Frobenius25,  one  can  show  that  in  physical 
variables  the  leading  behavior  is 


-(tf  -o 


2(11+1) 

(Sn-a(n+l) 


E  1 

^2  —  ^1  (tf  —  t) 


(64) 

(65) 


Here  |3  =  pj/p. 

Thus,  there  exist  two  solutions  (64-65)  for 
symmetric  perturbations  of  the  process  of 
explosive  boiling  of  a  liquid  droplet.  The  stability 
of  the  former  depends  on  the  ratio  of  the  densities 
of  the  vapor  and  the  host  fluid  a/p.  For  low 
density  of  the  host  fluid  (a>p)  the  solution  (64)  is 
also  unstable  and  grows  faster  than  (65).  When 
P>a  only  low  wave  numbers  are  unstable  in  (64). 
For  high  density  of  the  host  liquid,  the  solution 
(64)  is  stable.  That  means  that  the  droplet 
breakup  is  easier  in  the  gas  medium  than  in  the 
liquid,  which  is  a  physically  reasonable  result. 

On  the  other  hand,  the  solution  (65)  is 
always  unstable.  One  can  therefore  conclude  that 
the  process  of  the  explosive  boiling  of  a  liquid 
droplet  is  unstable. 


Conclusions 

In  this  paper  we  analyzed  the  stability  of 
explosive  boiling  of  a  liquid  droplet.  Two  related 
problems  were  considered.  We  showed  that 
evaporation  of  a  highly  superheated  droplet  is 
stable  and  we  studied  the  stability  of  a  thin 
expanding  liquid  shell.  The  results  of  the  former 
problem  can  be  also  used  for  a  contracting 
spherical  flame  while  investigation  of  the  latter 
may  throw  some  light  on  the  process  of  droplet 
breakup  at  the  final  stages  of  the  boiling. 

The  results  of  these  related  problems  were 
used  to  justify  the  assumptions  necessary  in 
studying  the  general  case  of  explosive  boiling 
stability.  It  was  found  that  the  process  is  unstable 
as  indeed  was  observed  in  the  experiments  of 
Shepherd  &  Sturtevant3.  The  instability  was 
obtained  for  all  wave  numbers,  which  again  is 
consistent  with  the  observation  of  Shepherd  & 
Sturtevant  (Ref.  3,  p.  388)  that  the  roughening  of 
the  bubble  surface  occurs  on  many  length  scales. 
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Fig.  1.  A  highly  superheated  evaporating  droplet. 
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Fig.  3.  Explosive  boiling  of  a  liquid  droplet. 


